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A new recursion procedure for deriving renormalized perturbation expansions for the one-dimensional anharmonic 
oscillator is offered. Based upon the ?i-expansions and suitable quantization conditions, the recursion formulae obtained 
have the same simple form both for ground and excited states and can be easily applied to any renormalization scheme. 
As an example, the renormalized expansions for the sextic anharmonic oscillator are considered. 

1. Introduction 

In the past few decades intensive investigations have been carried out on the one-dimensional anharmonic oscil- 
lator because of both its role in the modeling of quantum field theory and its usefulness in atomic, and molecular 
physics The conventional way to study the energy eigenvalues and eigenfunctions of this bound-state prob- 
lem is the practical application of perturbation expansions in a coupling constant. But the nonconvergence 
of obtained expansions ^'^ compels us to resort to modern procedures of summation of divergent series. Ones 
of the most common among them are various versions of the renormalization technique The inclusion of 
one or more free parameters corrects for the above mentioned deficiency, by controlling and accelerating the 
convergence of the expansions not only for eigenvalues but for eigenfunctions as well 

However, to provide a reasonable accuracy this method needs to know high orders of the perturbation series. 
Unfortunately we can easily obtain the high-order corrections only in the case of ground states when the simple 
recursion relations of the logarithmic perturbation theory i^-is -j-j^g Bender- Wu method ^ for performing the 
Rayleigh-Schrodinger perturbation expansion are applied. The description of radially excited states involves 
including the nodes of wave function in consideration that makes these approaches extremely cumbersome and, 
practically, inapplicable. 

Very recently, a new semiclassical technique for deriving results of logarithmic perturbation theory within 
the framework of the one-dimensional Schrodinger equation has been proposed Based upon an ?i-expansion 
this technique leads to recursion formulae having the same simple form both for the ground and excited states. 

The object of this paper is to extend the above mentioned formalism and to adapt it to the treatment of any 
renormalization scheme in terms of handy recursion relations within the framework of the united approach. To 
this end the paper is presented as follows. Section 2 contains a general discussion and the necessary assumptions 
$H in the semiclassical treatment of logarithmic perturbation theory. In Sect. 3 quantization conditions obtained 
are used for deriving the recursion formulae for perturbation expansions. Sect. 4 focuses on a derivation of the 
recursion formulae for renormalized perturbation expansions. Sect. 5 gives the example of explicit application 
of proposed technique. Finally, Sect. 6 contains concluding remarks. 



>< 



2. Semiclassical treatment of logarithmic perturbation theory 

In this section, our goal is to describe a straightforward semiclassical approach to obtaining results of the 
logarithmic perturbation theory. 

We shall confine our interest to bound-state problem for the Schrodinger equation in one degree of freedom 

- — U"ix) + V{x)Uix)=EU{x). (1) 
2m 

We shell study the discrete piece of the spectrum related to the potential function, V{x), having a single 
simple minimum and can be given by expression 

Vix) = ^mcj^x^ + J2f^^'^^- (2) 

i>i 

ipACS Numbers 03.65G, 03.65S 
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1 



However, if the potential has a finite number of local minima, the proposed below technique also can be 
applied. In this case the obtained eigenvalues are the energies of quasidiscrete states in a vicinity of each local 
minimum when the tunneling can be neglected. 

We want to restore the results of logarithmic perturbation theory for the nth eigenfunction and corresponding 
energy eigenvalue by means of the explicit scmiclassical expansions in powers of the Planck constant, h. 

Indeed, after the scale transformation, x Vhx, the coupling constants, fi, appear in common with powers 
of Planck's constant. Therefore, the perturbation series must be scmiclassical ?i-expansion, too. 

It has been proved that the energy eigenvalues under consideration should be concentrated near the 
minimum of the potential and should behave as 

E ^nLu{n+ 1/2) + J2EkiuJ,n)n^. (3) 

fc>2 

A few procedures are elaborated for computing the coefficients Ei{uj,n). They involve, in particular, applying 
the methods of the comparison equation '^^ and complex 'sprout' an analytic continuation in the ?i-plane 
various approaches within the framework of the WKB-approximation ^4-28. quantization using the methods of 
classical mechanics 26,29. g^j^j^ lastly, expansions in the ?i^^^-series 

Here we describe a new, simpler and more straightforward scmiclassical technique It involves two suc- 
cessive stages: first, we imply that the eigenvalue selection is effectively achieved by enforcing not the square 
integrability of the eigenfunction, but its analyticity, with applying the quantization condition; and second, we 
state such a classical analog of the system, in the limit ?i — > 0, which is appropriate for consideration of the 
low-lying levels. 

To begin with, we use the substitution, C{x) — hU' (x) /U{x), accepted in logarithmic perturbation theory, 
and rewrite Eq. (|l|) in the Riccati form 

hC'{x) + C'^ix) = 2m[V{x) - E]. (4) 

It will be recalled that we attempt to solve this equation in a scmiclassical manner with the following asymptotic 
series expansions 

CXD OC 

= c{x) = Y,Ck{x)h\ (5) 

fc=0 fc=0 

On substituting these expansions into the Riccati equation and collecting coefficients of the same powers of 
h, one obtains the system 

Cl{x) = 2mV{x), 
C'Q{x) + 2Cf){x)Ci{x) = -2mEi, 

k 

C'k_l{x)+Y,C^{x)Ck-^{x) - -2mEk. (6) 

1=0 

For nodeless states this system coincides with those derived by means of the standard perturbation approach 
and can be solved straightforwardly. However, when radial excitations are described within the framework of 
logarithmic perturbation theory the nodes of wave functions are included in some separate factor and considera- 
tion becomes extremely cumbersome. We intend to circumvent this difficulty by making use of the quantization 
condition, inspired by the WKB-approximation 3i,32_ 

Generally speaking, such the quantization condition simple means monitoring whether eigenfunction U{x) 
is analytic inside a closed contour. We would remind that a necessary condition for analyticity is 

27ri / U{x) ' ^ ' 

which must be a non-negative integer, giving the number of zeros of U(x) which lie inside the contour of 
integration. Hence, the transition to the logarithmic derivative, C{x), allows the nodes of wave functions to be 
easily involved in consideration. 
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Taking into account that due to the Sturm-Liouville theory the eigenfunction of the nth radially excited 
state has precisely n zeros on the real axis, we have 

(j) C{x)dx = ^^^h'' (j) Ck{x)dx = nn, n = 0,l,2, ... (8) 

k—O 

where a contour of integration encloses the nodes of wave function and no other singularities. 

The above quantization condition must be supplemented with the rule of achieving a classical limit for 
its right-hand side. Unfortunately, under the influence of standard textbooks on quantum mechanics, the 
semiclassical 7i-expansions are usually associated solely with the WKB-approach, for which the rule of achieving 
a classical limit is 

7i — + 0, n — > cxD, hn — const. (9) 

That implies that the WKB method is more suitable for obtaining energy eigenvalues in the limiting case of 
large quantum numbers. 

The semiclassical treatment of the low lying levels demands the use of the different, as compared with the 
WKB-approach, rule of achieving a classical limit, namely 

n-O(l), hn-^O. (10) 

This rule does provide lowering a particle to the bottom of a potential well as 7i — + and was formerly 
applied in deriving coefhcients of the 1/N-expansions 33-36^ ^^^i be shown below, in one dimension case 

this way lead us to restoring the logarithmic perturbation series. 

Thus, in view of the new rule of passing to the classical limit (|o|), the quantization conditions are 
rewritten in the final form as 

— (f Ci(x) dx ^ n, — (f Ck(x)dx^O, fc ^ 1. (11) 

2tTI J ZTTi J 

A further application of the theorem of residues to the explicit form of functions Ck{x) easily solves the 
problem of taking into account nodes of the wave functions for excited states. 

3. Recursion formulae for perturbation expansions 

We proceed now to deriving the recursion relations for obtaining the nth eigenfunction and corresponding energy 
eigenvalue. 

Let us consider the system (^) and investigate the behaviour of the functions Ck (x) . From the first equation 
it is seen that Co(x) can be written in the form 



Co{x) ^ -[2mV{x)]^^^ = -mux{l + ^ V /.x^'/' = 2; V C°x\ (12) 

i>l 1=0 

where the minus sign is chosen from boundary conditions, and coefficients C? are defined by parameters of the 
potential through the relations 



Cq = —mu!, 



2m(jj 

p=i 



From (H) we recognize that Co(0) = and, consequently, the function Ci{x) has a simple pole at the origin, 
while Ck{x) has a pole of the order of {2k — 1). Thus, Ck{x) can be represented by the Laurent series 

00 

Ck{x)=X^-^''J2'^^^'^ '^^l- (14) 
j=0 



With definition of residues, this expansion permits us to express the quantization conditions (11) explicitly in 
terms of the coefficients C,^ as 
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= nSi,k, (15) 

where the Kronecker delta is used. Thereby the common consideration of the ground and excited states has 
indeed proved to be possible. 

The substitution of the series (^ and (|l|) into (|) in the case i ^ 2k — 2 yields the recursion relation for 
obtaining the coefficients Cf : 

= -^[(3 -2k + ^)Ct' + E E ^i^'-P + 2E (16) 

j = l p=0 p=l 

whereas if i = 2fc — 2 we would find the recursion formula for the energy eigenvalues 

A; 2fe-2 
j=0 p=0 

It is easy to verify that proposed technique does indeed restore the results of the logarithmic perturbation 
theory. Upon calculating first corrections to the energy eigenvalues we arrive at the form familiar from standard 
textbooks 



E = h^{n + 1) - r^'^An' + + n'^,{n' + n + 1) + 0{h% (18) 



with the obvious constraints 

The equations (^6|) and ( [l7| ) have the same simple form both for ground and excited states and provide, in 
principle, the calculation of the perturbation corrections up to an arbitrary order in the analytical or numerical 
form. 

Another salient feature of this technique is the fact that for obtaining corrections to the energy eigenvalues 
and eigenfunctions we need not have initially the explicit form of the exact solution for the unperturbed potential. 
This solution is found automatically in the process of solving recursion formulae. 

Moreover, in this way the quasi-exact solutions may be obtained as well. Indeed, to be more concrete, let 
us consider the ground state for the anharmonic potential 

V{x) ^V2X^ + Vi + Ve (20) 

Then the lowest order corrections derived by means of formulae ( p^ and ( [l7| ) are written as 

El = v/^, 
3T4 



E2 — 

E-t = 



2V2' 

3 {-IVj^ + 5V2V6) 
3 



-9 {~37V4'' + A0V2ViVe) 
Ei — 1 , (21) 



E. 



■5 



E, 



6 



setting 



-15 (2059 Vi - 2992 V2 Va + 466 V-i ) 
641^2"^^ 

9 (101859 - 186380 V2 Vi Ve + 61420 V-i Vi V^,- 
128W^ 



^2 = yli{^ y^) - 3 (22) 
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immediately leads to the analytical expression for the ground state energy and the wave function ^"^ 

E = Vi/{2^Q) (23) 

and 

*(.t) - exp(-F4 xV(4/t^) - a;V4) (24) 
where for the simplicity we put fi = 2m = 1. 



4. Renormalized perturbation expansions 

The renormalization of perturbation expansions is based on a redistribution of values of series items for a partial 
summation of them and improving the convergence of expansions. From mathematical point of view it means the 
transition from the asymptotical expansion according to Poincare to the asymptotical expansion according to 
Erdelyi Characteristic of this method is that the zeroth approximation involves a set of artificial parameters, 
which are not contained in the original Hamiltonian and are determined order by order after calculating the 
physical quantity perturbatively. 

For the discussed one-dimensional anharmonic oscillator, as a trial parameters we can choose a frequency 
and a mass. Though proposed technique is easily applied to the mass renormalization we do consider only the 
case of renormalization of a frequency as being more physically motivated. For this purpose it is enough to 
think of the harmonic oscillator frequency incoming in the potential as a function of Plank's constant variable, 
with subsequent its expansion in an 7i-series However, for later use it is more convenient to take 



oo 

,2 _ 



(jj 

k=0 



The recursion system (||) then reads 



C^(.t) = m^ujlx^ + 2to E fi- 



i>l 

C'oix) + 2Ca{x)Ci{x) = ra^uj\x^ - 2mEx, 



i+2 



i=0 

with 



fc 

C'k^l{x) + Y,C^ix)Ck-^{x) = m'^ujlx'' -2mEk, (26) 



(27) 



j^k 

^ " 2Co" 



° p=i 

In the case A: > 1, when i ^ 2k — 2 and i ~2k — 2 that results in 

-m^uJl6^ak + {i-2k + t)C^-^ 

rv— J. t i 

E E '^p'^i-p + 2 E ^p^i-p 
j = l p=0 p=l 

k 2fe-2 

2mE, = ~C',-},~J2J2^P^2k^2-p, (28) 

j=Q p=0 

respectively. Here the coefficients are defined by the chosen version of renormalization. 

In practice, the one-parameter schemes are usually used. They are obtained with the truncation of series 
(H) as 

cj^ = oj^„ + u;ln'' , (29) 



fc-i 
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where ujq is a trial frequency and an order in h of the remainder is determined by the anharmonicity of a 
potential (||). 

A free parameter, loq, is order-depended. At any finite order of approximation, N, its optimal value is picked 
up by the condition of "minimal difference" En{luo) = 39. 4o^ g^j^j "minimal sensitivity", {d/dujo) J2f=o ^ii'^o) = 
^^''^^ , because perturbatively calculated quantity should not depend on the artificial parameter at all. 

Notice that our approach can be also applied to more general scheme such as, for instance, the renormgroup 
conditions '^^ or the simple equating to zero some corrections: 

E,{u;o,...,ujk)=0, i = 2,3,...,N. (30) 



5. Example of application 



As a rule, various schemes of renormalization are illustrated with describing ground states of the quartic anhar- 
monic oscillator. Here we consider the less traditional example of obtaining energy eigenvalues for both ground 
and excited states of the sextic anharmonic oscillator with a potential 



(31) 



where we put tu = m = 1 for simplicity. 

The analytical representation of perturbation expansion for its energy obtained by using formulae (|l^) and 
ITn has the form 



El 



1 



A (1 + 2n) (3 + 2n + 2n^) , 
-2"^ (1 + 2n) (3495 + 4538 n + 5324 + 1572 + 786 n^) , 
2-"5A3 (l + 2n) (247935 + 444014 n + 600050^2 + 323868 



-191424 + 35388 + 11796 n^) 



(32) 



This series diverges for any finite value of A and requires the use of some renormalization procedure. Since for 
discussed potential the above series is, in fact, the expansion in powers of the Planck constant squared, the trial 
frequency within any one-parameter scheme of renormalization is given by (ES) as 



(33) 



By taking use of relations (^8|), the first terms of the renormalizcd perturbation series for the energy eigen- 
values then read 



El = 

E3 - 
E5 = 



Er 



2-^luq^ {l + 2n) [5A {3 + 2n + 2n^) +Auj'^ (l-^^o)] 



-2"*W(7'^ {l + 2n) X 



X A^ (3495 + 4538 n + 5324 + 1572 + 786 n*) 
-Hl20Aw^ (1-wo^) {3 + 2n + 2n^) + 16 ujo^ {l - uja^ 



2~"wo" (l-|-2n) X 



5 A^ (247935 + 444014 n + 600050 + 323868 ) 



191424 + 35388 n'' + 11796 n'^ 



+28X^ujo^ (1 
X (3495 + 4538 n + 5324 + 1572 



786 71-* 



-t-1200Awo'' (1 -cjo^)^ (3 + 2n + 2n2) + 64, ujq^ (l - ujo^Y 



(34) 



where UJ2 is excluded due to equation (B3). 
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The obtained recursion formulae enable us to calculate high orders of the renormalized perturbation series 
for wide ranges of the coupling constant A and the radial quantum number, n. 

We have computed the partial sums of N corrections to the energy eigenvalues with defining the free 
parameter ujq, order by order under the conditions of the minimal difference, £'jv(a;o) = 0, and the minimal 
sensitivity, {d/d(jJo)E^{(jJo) — and {d/dojo)^^^^ Ei{uJo) = 0. 

Contrary to the current view we recognize that in the case of sextic anharmonic oscillator all these 
schemes appear to give approximately an equal accuracy. 

It has been observed that in the range of values < A < 1000 and < n < 1000 the best accuracy can 
be achieved in the limiting case of the weak coupling, for small A, n. The obtained energy eigenvalues are 
most inaccurate in the intermediate region (A = 0(1) and < n < 2). Further, the Nth-order renormalized 
approximation becomes improve as n and A increase. 

The typical behaviour of the spcc;d and accuracy of these one-parameter schemes is illustrated by Table 1. 
Here the partial sums of N renormalized corrections to energy eigenvalues derived under the condition of mini- 
mal sensitivity, {d/dcoo) X]t=o ^ii^o) — 0, are compared with the energy eigenvalues obtained by the numerical 
integration of the Schrodinger equation. As the numerical integration procedure it was used the improved shoot- 
ing method with the Noumerov finite difference scheme, described in literature which ensures the necessary 
exactness of calculation for the Sturm-Liouville problems. Notice that the condition of minimal sensitivity yields 
a large number of the trial wo-values from which we choose the one with the flattest extremum ^2,44 

6. Conclusion 

A new simple recursion procedure, suitable to applying any renormalization scheme of improving the per- 
turbation expansions for the bound-state problem of anharmonic oscillator within the framework of the one- 
dimensional Schrodinger equation, has been proposed. 

Based upon the scmiclassical treatment, including ?i-expansions and proper quantization conditions, the 
main disadvantage of logarithmic perturbation theory, namely, the cumbersome description of excited states 
when nodes of wave functions are taking into account, has been avoided. 

The derived formulae have been adapted for the use of any renormalization scheme. As a result, the recursion 
formulae, obtained for renormalized perturbation expansions, have the same simple form both for ground and 
excited states and provide, in principle, the calculation of renormalization corrections up to an arbitrary order. 

Moreover, the advantage of the method is such flexibility that within the united approach one can apply 
various schemes of renormalization and choose different values of the trial parameter depending on the order of 
perturbation. 

The proposed technique has been developed for investigation of the bound-state problem within the one- 
dimensional Schrodinger equation with any potential described by an analytic function having a simple global 
minimum. However, taking into account that the potential having a finite number of local minima becomes, as 
?i ^ 0, a sum of infinite deep wells, our procedure may be applyed as well for investigation of bound states in 
a vicinity of a local minimum, when tunneling can be neglected. 
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Table 1; The sequences of the partial sums of N renormalized corrections to the energy eigenvalues of the sextic 
anharmonic oscillator with the potential V{x) = | [x^ + Aa;^), computed within the one-parameter scheme of 

rcnormalization under the condition of minimal sensitivity, (d/dujo) X^^q EiiuJo) = 0. The values of Enum were 
obtained by numerical integration of the Schrodinger equation. 



N 


n = 


n= 1 


n = 5 


A = 0.01 


A = 10 


A = 0.01 


A = 10 


A = 0.01 


A = 10 


1 


0.508693705 


1.161458 


1.55611747 


4.210051 


6.59434725 


25.95659 


3 


0.508378396 


1.110292 


1.55399477 


4.054344 


6.57502024 


25.54466 


5 


0.508371342 


1.104354 


1.55397174 


4.047270 


6.61788050 


26.44265 


10 


0.508370692 


1.102706 


1.55398991 


4.056856 


6.61763448 


26.42384 


15 


0.508370674 


1.102541 


1.55398998 


4.057586 


6.61764001 


26.42596 


20 


0.508370673 


1.102651 


1.55398999 


4.057838 


6.61764261 


26.42806 


25 


0.508370689 


1.102586 


1.55398995 


4.057637 


6.61763908 


26.42459 


30 


0.508370674 


1.102729 


1.55398996 


4.057495 


6.61763913 


26.42449 


35 


0.508370726 


1.102819 


1.55398997 


4.057749 


6.61764005 


26.42450 


40 


0.508370676 


1.102768 


1.55398996 


4.057442 


6.61763907 


26.42474 


45 


0.508370682 


1.102796 


1.55398996 


4.057401 


6.61763907 


26.42484 


50 


0.508370681 


1.102829 


1.55398996 


4.057410 


6.61763906 


26.42479 
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